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A generalized version of the valence-bond Monte Carlo method is used to study ground state properties of the 1 + 1 dimensional 
quantum (2-state Potts models. For appropriate values of Q these models can be used to describe interacting chains of non-Abelian 
anyons — quasiparticle excitations of certain exotic fractional quantum Hall states. 
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1. Introduction 

The notion of a valence bond, a simple singlet state formed 
by two localized spin- 1/2 particles, captures much of the 
physics of how pairs of electrons are correlated in the real 
world. This is certainly true when describing chemical bonds 
in molecules, where a local picture of electron correlations is 
clearly appropriate; however, the language of valence bonds has 
also proven useful for describing a variety of possible singlet 
ground states of quantum spin systems in the thermodynamic 
limit, including resonating valence-bond states and valence- 
bond solids. 

Valence bonds can also be used for Monte Carlo simula- 
tions, as first shown in the context of variational Monte Carlo by 
Liang, Doucot and Anderson in 1988 [1 1. More recently, Sand- 
vik (2J has introduced a projector Monte Carlo method known 
as valence-bond Monte Carlo ( VBMC) which can be used to ef- 
ficiently sample ground states of quantum spin systems directly 
from the valence-bond basis. 

In this Proceedings we describe some details of the meth- 
ods used in our recent work on VBMC simulations of both uni- 
form and random spin- 1/2 antiferromagnetic Heisenberg chains 
and the closely related 1 + 1 dimensional quantum Q-state Potts 
models 0. One motivation for this work is that for certain 
values of Q these models describe interacting chains of non- 
Abelian anyons |0], exotic quasiparticle excitations which are 
thought to arise in certain fractional quantum Hall states 0. 
The results presented here are all for the case of uniform mod- 
els (i.e., with no disorder) for which a number of exact results 
are known which can be used to benchmark the method. 



2. Hilbert space and models 

We begin by describing the valence-bond basis for a chain 
of N spin- 1/2 particles. This basis is made up of valence- 
bond states — states in which all N particles are paired up to 
form N/2 valence bonds. Figure [TJ a) shows two (normalized) 
valence-bond states \a) and \/3). Both these states are exam- 
ples of non-crossing valence-bond states, meaning that no two 
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Figure 1: (a) Two normalized non-crossing valence-bond states \a) and |/3) in 
which pairs of particles (solid circles) are connected by valence bonds (solid 
lines), (b) To determine the overlap of these states one overlays the two valence- 
bond configurations and counts the number of closed loops, Ni oops . The overlap 
is then (a\f3) = d Nh "P s ^ N ^ 2 where d = 2 for the case of spin- 1/2 particles. 



valence-bonds cross each other, or, equivalently, that the total 
spin of all the particles between any two particles connected by 
a valence bond must be 0. 

For ordinary spin- 1/2 particles the set of non-crossing 
valence-bond states forms a complete and linearly independent 
basis spanning the space of all total spin states [6|. The num- 
ber of these non-crossing states, and hence the dimensionality 
of the total spin Hilbert space, grows asymptotically as 2 N for 
large N, as one would naturally expect for N spin- 1/2 particles. 

One price to be paid for doing numerical calculations with 
the valence-bond basis is that it is a nonorthogonal basis. The 
rule for determining the overlap of any two valence-bond states 
is shown in Fig. [TJb) — one simply overlays the two valence- 
bond configurations and counts the number of loops formed, 
Nioops- The overlap is then (a\/3) = d N,mp, ~ N ' 2 , where, for spin- 
1/2 particles, d — 2. 

When d + 2 the Hilbert space spanned by non-crossing 
valence-bond states is still perfectly well defined provided 
d > 1, although it no longer describes a system of ordinary 
spin-1/2 particles. For certain values of d, specifically when 
d — 2 cos ^ where A; is a positive integer, this Hilbert space can 
be interpreted physically as describing the 'topological charge 
0' sector of a system of N non-Abelian anyons described by 
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Figure 2: (a) Diagrammatic representation of a singlet projection operator, and 
(b) action of two projection operators on a given valence-bond state. In the first 
case (top), the loop formed introduces a factor of d leading to an overall factor 
I = 

factor is \ . 



su(2)k Chern-Simons-Witten theory Q. For these special val- 
ues of d, the non-crossing valence-bond states are no longer lin- 
early independent and the dimensionality of the Hilbert space 
for N particles can be shown to grow asymptotically not as 2 N 
but as d N . The quantity d is known as the quantum dimension 
of the particles Q. 

Such su(2)k anyons are thought to arise physically in cer- 
tain experimentally observed fractional quantum Hall states, 
notably the state with Landau-level filling fraction v = 5/2 (cor- 
responding to k — 2) and, possibly, the v = 12/5 state (corre- 
sponding to k — 3) |5|. 

Having defined the relevant Hilbert spaces we now turn to 
the model Hamiltonians studied here. To define these models 
we first describe the action of a nearest-neighbor singlet pro- 
jection operator in the valence-bond basis. Let n° denote the 
singlet projection operator acting on sites i and Figure 
|2ja) shows a useful diagrammatic representation of this oper- 
ator. Figure |2|b) uses this representation to illustrate the two 
distinct cases which can occur when acting on a valence-bond 
state with n°. Either one applies the projection operator to two 
sites which are connected by a valence bond, in which case the 
projection operator has no effect on the state, or one acts on 
two sites which are each connected to different sites by valence 
bonds, in which case the projection operator forms a singlet be- 
tween the two particles it acts on, as well as the two particles 
connected to them, and gives an overall factor of 1 Id. 

The Hamiltonians we consider here all have the form 



(1) 



For ordinary spin-1/2 particles with d — 2, the singlet projec- 
tion operator can be expressed as n° = j - S, • S,+i and cor- 
responds to an antiferromagnetic nearest-neighbor Heisenberg 
chain. More generally, for arbitrary d, one can readily check 
that the operators Ui = d n° satisfy the so-called Temperley- 
Lieb algebra |8|, 



Uf = dUi, 
UiU i±l Ui = u it 
[Ui,Uj] = 0, \i-j\>l. 



(2) 



This algebra appears in the study of the 2-dimensional Q-state 
Potts model with Q = d 2 , and as a consequence it can be shown 



that the models ([T]i are equivalent to the 1 + 1 dimensional quan- 
tum (2-state Potts models [9]. Furthermore, for the special val- 
ues d = 2 cos these models correspond to a sequence of con- 
formally invariant Andrews-Baxter-Forrester iflOl models with 
central charges q = 1 - 6/(k + l)(k + 2) ifTTI . As stated above, 
these models can be thought of as describing chains of interact- 
ing non-Abelian anyons 0). 



3. Valence-bond Monte Carlo 

The basic idea behind VBMC is to act on a particular 
valence-bond state \S) repeatedly with —H in order to project 
out the ground state. The results of this projection after n itera- 
tions can be expressed as follows, 



{-Hf\S)= Yj U l '•' n °J S >- 



(3) 



Since acting on a non-crossing valence-bond state yields an- 
other non-crossing valence-bond state, (see Fig. ^b)), it fol- 
lows that 
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where \a) is a non-crossing valence-bond state with the same 
norm as \S) and /I;, = d~ m where m is the number of time a 
projection operator acts on two sites which are not connected by 
a valence bond in the process of projecting \S ) onto the state \a). 
Figure [3] shows diagrammatic representations for two different 
terms of the form Q with N - 12 and « = 6 (in our simulations 
we study system sizes up to N = 1024 and take n = 2QN). 

In the end the projection ([3]) leads to an expression for the 
ground state \if/) of the form 



|<A> = J] w(a)\a). 



(5) 



In VBMC one samples the valence-bond states \a) contributing 
to with probability w(a)/ 2a w(JS) by updating the sequence 
of projection operators n? , • ■ ■ , IT? using the usual Metropolis 
method. One such Monte Carlo update is shown in Fig. [3] 

Given any observable O with expectation values 0(a) = 
(a\0\a) in the non-crossing valence-bond states |or), VBMC can 
be used to compute the average (O) = £ a w{a)0(a)/ £ a w(a) 
for any state \iff) of the form (|5j, provided w(a) > (which is the 
case here). In what follows, angle brackets will always denote 
this average, though it should be noted that (O) will in general 
not be equal to the true expectation value {if/\0\if/)/{i//\if/), both 
because the valence-bond states are nonorthogonal and because 
the weight factors w(a) are amplitudes and not probabilities. 
Of course the true quantum expectation value of any operator 
can be computed by VBMC if one carries out the projection on 
both the bra state and the ket state [2 |. However, here we focus 
on those quantities which can be calculated efficiently using the 
"one-sided" VBMC described above. 
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Figure 3: Diagrammatic representations of the action of two sequences of singlet projection operators on a given starting valence-bond state. In each diagram the 
resulting valence-bond state \a) is determined by the open loops which terminate at the top of the diagram and the amplitude Xi w ..j s is the product of the factors 
listed on the right of the diagram. Each of these factors is either 1 or 1 /d, depending on whether the corresponding projection operator, listed on the left of the 
diagram, forms a closed loop or not. The two diagrams shown can be viewed as 'before' and 'after' pictures for a single Monte Carlo update in VBMC in which 
the third projection operator from the bottom is shifted from rtj| to n[J and the update is accepted or rejected according to the usual Metropolis rule. The quantity kl 
which is used to compute the valence-bond entanglement entropy (see text) is defined to be the total number of bond leaving a given block of L sites (red rectangle). 



4. Results 



4.1. Ground state energy 



One quantity which is easily computed using one-sided 
VBMC is the ground state energy Eq. The procedure given in 
|2| for calculating for spin- 1/2 systems can be trivially gen- 
eralized for arbitrary d and leads to the following expression, 
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Here, w, is equal to 1 if sites i and i + 1 are connected by a 
valence bond and l/d if they are not for a given valence-bond 
state |or), and, as described above, the angle brackets denote an 
average over these valence-bond states weighted by the ampli- 
tudes w(a). 

In fact, the ground state energies of the models ([T]i can be 
found exactly. This can be seen by noting that for any d the 
Temperley-Lieb operators can be represented using spin- 1/2 
operators as U t = 2(SfS* +1 + S y S y +l ) + d(l/4 - + 
i Vl -<i 2 /4(Sr +1 - SI) |9). For the case of open boundary con- 
ditions the models ([TJ can then be mapped onto spin- 1/2 XXZ 
chains with external (non-Hermitian) fields applied to the two 
ends (the staggered field term in the expression for U, cancels 
in the "bulk" of the chain). In the thermodynamic limit, the 
ground state energies will not depend on boundary conditions, 
and the values of Eq for the models ([T]i with periodic bound- 
ary conditions should be the same as that for the corresponding 
XXZ models. Using the expression for the ground state energies 
of the XXZ models found using Bethe ansatz by Yang and Yang 
Ifl2l it is straightforward to obtain the following expression for 
the ground energies of the models ([TJ, 
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Figure 4: Ground state energy per site Eq/N as a function of quantum dimen- 
sion d. The solid red line is the exact Bethe ansatz result, and the blue squares 
are the results of our VBMC simulations (error bars are smaller than symbol 
size). 



Figure |4] shows the ground state energies we obtained by 
evaluating the expression (|6jl by VBMC for d = 2 cos ^ with 
k = 2, 3, 4, 5, 6 and oo. The red line is the exact energy (|7]l as a 
function of d. The fact that our numerical results clearly agree 
with the exact Bethe ansatz results should be seen as evidence 
that VBMC can indeed be used to simulate the models ([TJ with 
arbitrary d. 

4.2. Bond length distribution and valence-bond entanglement 

Another quantity which is natural to compute using one- 
sided VBMC is the so-called valence-bond entanglement en- 
tropy Ifl3l [T4l . If « L is defined to be the total number of bonds 
leaving a contiguous block of L sites in a given valence-bond 
state \a) (see Fig. [3]) then the valence-bond entanglement en- 
tropy is defined to be S (L) = (ni). While this quantity was orig- 
inally only defined for spin- 1/2 systems, the definition clearly 
generalizes to the models considered here for arbitrary values 
of d, as first noted by Jacobsen and Saleur [15 |. In this same 
paper, Jacobsen and Saleur also obtained analytic results for the 
L » 1 scaling of S (L). For all d < 2 they found that S (L) scales 




Figure 5: Semi-log plots of the valence-bond entanglement entropy S(L) as functions of the block size L c , and log-log plots of the bond-length distribution P(l) 
as a function of bond length l c for d = y2 and d = 2. Here L c = (N/n) sm(nL/N) and l c = (N/n) sm(jrl/N) are the so-called conformal lengths which we use to 
minimize finite size effects when L and / are nearly equal to half the system size. The solid lines in these plots are the analytic results for the asymptotic scaling of 
S(L) and P(l) which follow from the results obtained in 1 15 1. Results are for chains with N = 1024 sites and the parameter n is taken to be 20N. 



logarithmically with L with a d dependent coefficient. It should 
be noted that this exact result was based on a mapping to a long- 
wavelength field theoretic description of the models, and so it 
is worthwhile testing this prediction numerically. 

Before describing our results for S(L), we note that this 
quantity is closely related to the bond-length distribution, P(l), 
for valence bonds. This distribution is defined to be f(Z) = 
(bi)/(N/2) where bi is the number of valence bonds of length / 
in a given valence-bond state \a). It is readily shown [16] that 
P(l) is related to S (L) by, 

N 

S(L) = 2YP(l)tsan(l,L). (8) 
1=1 

(Note that P(l) = for even I.) This expression, together with 
the result of lfl5ll that S(L) scales logarithmically with L for 
d < 2 implies that, for these values of d, P(l) should follow an 
inverse-square power law for I » 1. 

Figure |5] shows our VBMC results for S(L) and P(l) for the 
cases d — V2 (corresponding to the critical one-dimensional 
transverse field Ising model) and d = 2 (corresponding to 
the spin- 1/2 antiferromagnetic Heisenberg chain). The log-log 
plots of P(l) vs. / clearly demonstrate the predicted inverse- 
square power law dependence for the bond-length distribution. 
The semi-log plots of S (L) also show the expected logarithmic 
scaling of the valence-bond entanglement entropy, for L » 1 . 

5. Conclusions 

To summarize, we have shown that the VBMC method of 
Sandvik [2] can be straightforwardly generalized to study the 
1 + 1 dimensional quantum Q(= c/ 2 )-state Potts models. For 
d — 2 cos these models describe chains of interacting non- 
Abelian anyons, exotic quasiparticle excitations believed to ex- 
ist in certain experimentally observed fractional quantum Hall 



states |Hl|5l. The ground state energies, bond-length distribu- 
tions and valence-bond entanglement entropies of these models 
were computed using VBMC and compared to various known 
exact results. This work sets the stage for our recent VBMC 
study of the effect of disorder on these models (3). 
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